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We discuss the relevance of the Stillinger and Weber approach to the glass transition investigating 
the non-equilibrium behavior of models with non-trivial dynamics, but with simple equilibrium 
properties. We consider a family of ID constrained kinetic models, which interpolates between the 
asymmetric chain introduced by Jackie and Eisinger [Z. Phys. B84, 115 (1991)] and the symmetric 
chain introduced by Fredrickson and Andersen [Phys. Rev. Lett 53, 1244 (1984)], and the ID version 
of the Backgammon model [Phys. Rev. Lett. 75, 1190 (1995)]. We show that the configurational 
entropy obtained from the inherent structures is the same for all models irrespective of their different 
microscopic dynamics. We present a detailed study of the coarsening behavior of these models, 
including the relation between fluctuations and response. Our results suggest that any approach to 
the glass transition inspired by mean-field ideas and resting on the definition of a configurational 
entropy must rely on the absence of any growing characteristic coarsening pattern. 

PACS numbers: 75.10.Nr, 05.50.+q, 75.40. Gb, 75.40.Mg 



I. INTRODUCTION 

There is an old debate concerning the correct description of dynamics in the glassy stateR According to the general 
wisdom, undercooled liquids are in a locally equilibrated metastable phase, but fall completely out of equilibrium when 
the relaxation time exceeds the observation time. In this situation the glass ages and slowly relaxes towards equilib- 
rium. While it is widely accepted that the glass transition observed in laboratory is a purely kinetic phenomenon, it 
is still not clear whether a true (or what kind of) ergodicity breaking underlies the glassy behavior and whether the 
properties associated to a possible equilibrium transition manifest themselves on the experimentally accessible time 
scales. A particularly interesting problem concerns the precise mechanism leading to a slow relaxation and its relation 
with the ground states structure of the system. 

A possible description of the non-equilibrium regime is in terms of coarsening. The coarsening process is described 
by a length scale which grows in time driving the system towards equilibrium. The most typical scenario for a 
coarsening dynamics is found in a ferromagnet quenched down to a temperature below its critical temperature T c . 
After quenching, domains of positive and negative magnetization grow with time. The system acquires a macroscopic 
magnetization only when the typical domain size becomes of the order of the system size, leading to a nucleation 
process which involves overturning of large domains in favor of the dominant phase. Although it cannot be excluded 
that some type of coarsening behavior similar to that found in ferromagnets takes place in real glasses, there is no 
strong evidence, up to now, that any type of coarsening process occurs in the relaxation of an undercooled liquid. 

Another possible description for the observed non-equilibrium behavior calls for activated dynamics. The activated 
dynamics scenario is rather different from coarsening. No typical growing length scales are now present. The system 
approaches a disordered state, which has no correlation with the crystal state, via thermally activated jumps among 
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different configurations corresponding to structural rearrangements of spatially localized regions. In this scenario, the 
ordered crystal state has no special relevance. It occurs when fluctuations nucleate a crystalline droplet of size bigger 
than a given critical size, strongly dependent on external parameters such as temperature. Consequently crystallization 
can be completely inhibited by going to low enough temperatures. It is generally assumed that crystallization plays 
a role only for time scales much larger than those relevant for the relaxation of thermodynamic quantities in the 
undercooled phase (such as the enthalpy or specific heat). Indeed a glass transition also emerges for purely disordered 
systems where a crystal state does not existBB. This does not exclude the presence of heterogeneities in the glassy 
phase, whose role and influence on the dynamics still needs to be understoodO'El. 

Despite the enormous effort devoted to this subject it is still not clear how these two scenarios combine together 
into a final description of the glassy state. Moreover, if coarsening is the driving process for the relaxation of 
undercooled liquids, due to the completely unknown microscopic structure of the glass state (on which the system 
should asymptotically relax) it is unclear what should be experimentally measured in order to identify the growing 
domains. A similar problem is encountered in spin glassesD. Since a disordered state has not a periodic structure, 
a priori one does not know how to observe domains and the question of the nature of the domains in spin glasses is 
still unresolved!! However, for spin glasses such freezing of temporal spin pcarrelations leads to a divergence of the 
spin-glass susceptibility. Despite some experimentaB and simulation results&u, strong evidence for such a divergence 
is not found in structural glasses. The greatest difficulty in elucidating this issue is that a general non-equilibrium 
theory to deal with this class of systems is still missing and approximations to this problem remain partial. They 
usually— svork either in a limited range of time scales or in a limited range of temperatures (for instance, mode-coupling 
theory^). 

During the eighties a novel approach to the glass transition was proposed-by Stillinger and Webeio. This description 
of the undercooled liquid phase, inspired by the Adams-Gibbs-Di MarzicO theory, incorporates those features of the 
energy landscape relevant to the activated regime. The Stillinger- Weber approach is based on a decomposition of 
the configurational space into basins (also called inherent structures, hereafter referred to as IS) on the basis of the 
topology of the potential energy surface. This construction yields a definition of a configurational entropy in terms 
of dynamically accessible basi»fcand it is close in spirit to the equilibrium configurational entropy or complexity of 
mean-field spin-glass modelscJ'Ea. 

In this paper we study the relevance of the IS analysis proposed by Stillinger and Weber for non-equilibrium 
dynamics using one-dimensional constrained kinetic models. Kinetically constrained models were first proposed by 
Fredrickson and Anderscnllj in the attempt to provide a simple microscopic mechanism for understanding the purely 
dynamical transition predicted by the mode-coupling theory. In these models the slowing down of the dynamics is 
obtained through the introduction of dynamical constraints, compatible with detailed balance and Boltzmann-Gibbs 
distribution. What makes highly peculiar the relaxational behavior of these models is the fact that the slowing down 
of dynamics is only due to kinetic constraints, which prevent certain transitions from occurring. For a review on early 
results on these models see Ref. [l(| 

We find that the non-equilibrium properties of this class of models cannot be described in terms of the Stillinger aod. 
Weber approach, in contrast to what is found in otheunodels of structural glasses such as Lennard- Jones glassesoEZrt£l 
or finite-size fully connected disordered spin glassesE3. The main reason is that even though the dynamics is different, 
the configurational entropy derived from the Stillinger and Weber decomposition is the same for all models. While this 
conclusion is probably valid for all type of coarsening models, it is an open problem whether the precise mechanism 
for which the inherent structure entropy happens to coincide is more general and independent of coarsening behavior, 
i.e. it could also hold for higher dimensional constrained models with a different dynamics. 

The paper is organized as follows. In Section II we introduce the family of one-dimensional models studied. In 
Section III we discuss the Stillinger and Weber approach and the configurational entropy for these models. In 
Section IV we analyze and compare the non-equilibrium dynamics for the different models, in particular we consider 
the coarsening phenomena and the fluctuations response relation. Finally in Section V we present conclusions and 
discussion. Some technical points are analyzed in two Appendices. 



II. THE MODELS 



We consider two different classes of ID constrained models: the constrained Ising chain and the Backgammon (BG) 
model. For the former we shall mainly consider the two extreme cases of the symmetrically constrained chain (SCIC) 
introduced by Fredrickson and Anderser£j and the asymmetrically constrained Ising chain (ACIC) introduced by 
Jackie and Eisingeinj. The BG modelE3 is not a purely constrained kinetic model since there are not local constraint 
on the microscopic dynamics. The constraint here follows from the requirement that the number of particles remains 
constant. This induces a global dynamical constraint which slows down the dynamics as temperature is lowered. The 
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study of the BG model complement our investigation comparing its behavior with the locally constrained Ising chain 
models. 



1. The constrained Ising chain 
The Hamiltonian of the model is defined by 

v 

E = -£><, C 1 ) 

where <7j are Ising-like variables, which take the values 0, 1, and the index i runs over the sites of a D-dimensional 
lattice of volume V = L D . The model corresponds to a paramagnet in a field. The dynamics is of the Glauber type 
where the spins are randomly updated according to the following rule 

W{a, -» 1 - at) = 

with a and b = 1 — a positive real numbers. In this paper we consider the ID case D — 1 although the model is 
also interesting for larger dimensions. With exception of the configuration with all spins equal to 1 it is known that 
the space of configurations is an irreducible Markov chain, so that ergodicity is guaranteed and detailed balance is 
fulfilled. 

Depending on the values of a we may have different cases. In particular for a = b = | the model corresponds to 
the SCIC0, 

W(a l ^l-a t ) = ^(l-a l+1 +a l - 1 ) x rnm[l,exp(-/3A#)] (3) 

so that a spin can flip iff at least one of its neighbors is down. .— . 
If a is equal to or 1 then the model corresponds to the ACIC (a = 0)E2l, 

W(<Ti -»l-ffi) = mhi[l,cxp(-/3A£;)] x 5 ai _ lfi . (4) 

In this case a spin can flip iff its left neighbor is down, and the dynamics is more constrained than in the SCIC. For 
generic values of a the flipping of the spins may take a different probability depending if the left or right spin is down. 

The vast majority of works appeared in the literature focused on the previous two cases (0) and In the present 
paper, for completeness and for the scope of our discussion, we shall discuss when possible the behavior for generic a. 

The dynamics of these models cannot be solved exactly, even if several important r esult s are known, (i) The corre- 
lation time diverges Withe low temperature limit as t ~ exp(2/3) for the SCIC modelnJO and as r ~ exp(/3 2 / log(2)) 
for the ACIC modeloH (ii) In the SCIC model the.-S|taggered correlation functions relax exponentially fast with 
time and one can define two characteristic time scalesEH: the first one r = 1/(1 + exp(— /?)) non-activated, and the 
second one r\ — exp(/3) actuated. In particular the later defines the time scale before which no aging effects are seen 
in the coEcelation functionscll. (iii) In the SCIC model the hierarchy of dynamical equations can be exactly closed 
for T = OBEi In the Appendix A we show that this result, originally obtained for the SCIC model, can be easily 
extended to generic values of o. 



1 

D 



E( 

t=l,D 



a a 



i+fi 



bar-. 



min[l,exp(-/3A J B)], 



(2) 



2. The Backgammon (BG) model 
The energy (Hamiltonian) of the Backgammon model isEl 



A" 



E = -J2$n>,o, (5) 



i=l 



where n.$ = 0, 1, . . . , N counts the particles in each site (box) of a D-dimensional lattice of N = L D sites. The energy 
is given by the number of empty boxes. As before we consider the ID case, where D = 1. The dynamics is of the 
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Kawasaki type where the particles are randomly moved from one box to another and the change is accepted with 
probability 

W = min[l,exp(-/3AE)]. (6) 

Strictly speaking the BG model is not a constrained kinetic model since there are not local constraints on particles (or 
boxes). Nevertheless the conservation of particles number introduces a global constraint which makes the dynamics 
of the BG model glassy at low temperatures where a lot of particles accumulate in a small number of boxes and the 
further emptying of boxes becomes slower and slower as time goes on. 

In the original definition the dynamics of this model was meaa-field in the sense that particles could move from 
one box to any otheti-hax. This dynamics can be closed exactlyESI and many results have been obtained on its non- 
equilibrium behavioiESEJ. Here we are interested in the equivalent one-dimensional case where boxes are located on a 
closed ring (a chain with periodic boundary conditions) and particles can only move from one box to its left or right 
neighboring boxes. In this case the dynamics at low temperatures is driven by the coarsening of clusters of empty 
boxes similarly to that of spin-1 domains in the kinetically constrained Ising chain. Here, however, in addition to 
temperature-activated processes induced by energy jumps, the_system has entropic barriers which contribute in a non 
trivial way to the coarsening dynamics. In the original modelE3 with mean-field dynamics coarsening was absent and 
the entropic barriers were the only responsible for glassy dynamics leading in that case to an activated behavior. 

All these models share the common fact that, despite their dynamics, the thermodynamics is trivial and there are 
no equilibrium spatial correlations at any temperature. In particular, they do not show any finite-temperature phase 
transition. Therefore the non-trivial behavior follows only from the dynamics, which, dynamically constrained in the 
first case and ruled by entropic barriers in the second case, turns out to be glassy. 

In the next Sections we present a detailed investigation of the non-equilibrium dynamics of these models and discuss 
how their dynamics cannot be efficiently described in terms of an IS based configurational entropy approach. 

III. THE STILLINGER AND WEBER ENTROPY 
A. The Stillinger and Weber approach 

An interesting approach to investigate activated behavior in glasses was suggested in the eighties by Stillinger and 
Weberlill. This is based on the (natural) decomposition of the motion near the glass transition into intra-basin motion 
(within a valley) and inter-basin motion (jumps between valleys). In a "cage" picture the first motion corresponds 
to the motion of particles within a cage, while the second one to the creation or destruction of cages. This approach 
implements in a practical way the old idea that in the undercooled liquid a strong separation of timescales of the two 
motions occurs near the glass transition. The approach bears some resemblance to the Edwards packing entropy in 
the context of granular systemsE!rE3. 

Within the Stillinger and Weber (SW) approach each configuration of the system is mapped into a local minimum of 
the energy through a local potential energy minimization which start from the given configuration. The local minimum 
was called inherent structure (IS), while the set of configurations flowing into it defines the basin of attraction or 
valley of the IS. Following SW one then constructs a IS-based thermodynamics decomposing the partition function 
sum into a sum over IS with the same energyO 

Z(T)^"£p IS (e,T), (7) 

e 

with 

P IS {e, T) = exp N [-0e + s c (e) - 0f(0, e)] , (8) 

where s c (e), defined as the configurational entropy, yields the number of different IS with energy e: 0(e) = 
exp(7Vs c (e)). The term f(fi,e) accounts for the free energy of the IS-basin of energy e, i.e., the partition sum 
restricted to the basin of attraction of IS with energy e. In each IS-basin the energy has been shifted, so that the IS 
has zero energy, and / accounts only for energy differences. In general f(/3, e) may have a non-trivial dependence on 
the energy if the IS-basin of IS with different energy are quite different. When the temperature is such that only the 
states near the bottom of the IS-basin do contribute then it is reasonable to expect that /(/3, e) is roughly independent 
of alatj. Another case where the dependence of /(/?, e) on e is expected to be negligible is when thf IS-basins are 
narrow and contain few configurations. This approximation works very well for REM-likc modclsE3u3. When the 
e-dependence of / can be neglected the configurational entropy s c (e) can be obtained directly from (^). 
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As long as the configurations counted on the r.h.s. of eq. (|7h are the most relevant for the thermodynamics at 
temperature T the above construction is totally legitimate as far as thermodynamics is concerned. It is only a different 
way of summing the partition function. Nevertheless the s c (e) obtained with the SW decomposition is in some sense 
a dynamical quantity since the projection between configurations and IS basins can be seen as the zero-temperature 
dynamics of the model. For this reason we will refer to it as Stillinger- Weber configurational entropy to distinguish 
it from other possible definitions of the configurational entropy taken from mean-field conceptsoEZl. This poses the 
question, first raised by Monasson and Birolio, on the relevance of s c (e) and IS in general for dynamics. 

It is clear that once the energy and the rules of the dynamics are given the IS can always be defined. For example 
for spin-glass models with quenched disordered variables taken from a continuous distribution, the dynamics usually 
consists of Monte Carlo updates (for instance, single spin- flips). The IS are then identified as the final configurations 
reached after a sequence of Monte Carlo moves where the spin which yields the largest decrease of energy is identified 
and flipped. Consequently, IS are stable against single spin-flips but not for higher-number of spin flips. Biroli and 
MonassonE-J conclude then that IS are ill-defined because s c (e) depends on the number of spin-flips, which make the 
configuration stable (at least, for not- fully connected models). We disagree with this conclusions because, as noted 
above, the IS and the SW entropy are intimately related to dynamics, and therefore it is not a surprise that changing 
the dynamics the IS and the SW entropy in general change. What, in our opinion, is ill-defined is to speak of IS 
without specifying the dynamics. 

Nevertheless the question posed by Biroli and Monasson is far from being trivial. Indeed despite the fact that the IS 
and the SW are dynamical quantities it is far from obvious that they contain all (or almost all) relevant informations 
on the dynamics on long time scales. This is a well known problem in the theory of dynamical systems. The SW 
decomposition can be seen as a mapping of the true dynamics at a given temperature onto a symbolic dynamics given 
by the dynamics of the IS. The obtained symbolic dynamics gives-a good description of the original one only if the 
mapping defines what is called a "generating partition" , see e.g.E3. In general for a generic dynamics it is not at 
all trivial to demonstrate that such a partition exists, and even if it does exist, how to find it. We can then recast 
the question posed by Biroli and Monasson in the following way: does the SW decomposition lead to a generating 
partition, or at least to a good approximation of it, for the long time dynamics of glasses near the glass transition? 

In general we can answer to this question only a posteriori. We define a possible partition and then check if this 
reproduces the desired features of the dynamics. However we can try to find under which conditions the answer 
could be affirmative. Usually to find a generating partition, or a good approximation of it, a good starting point is 
by looking at the "physical" properties of the dynamics. The SW mapping replaces each configuration in a IS-basin 
with the IS itself. Therefore it is clear that this mapping will be a good mapping if the systems spends a lot of time 
inside the basin. Under this assumption the dynamics on time scales larger than the typical residence time inside a 
IS-basin should be quite well described by the IS dynamics. This scenario is typical of a many valley dynamics with 
activated dynamics. It is also clear from the above discussion that if the IS mapping is a good mapping it does not 
matter which configuration inside the IS-basin is used to represent the IS-basin. It may be the IS itself or any other 
configuration in the basin. In a recent study of finite-size mean-iield spin-glass models, which share the properties of 
structural glasses, this independence has been indeed observedEI 

On the contrary for dynamical processes described by a coarsening process this description should in general 
fail because dynamics proceeds through geometrically correlated configurations. Barrat, Burioni and Mezard have 
shownEJ that the difference between the two scenarios has a simple manifestation in how dynamical trajectories 
departing from the same configurations separate in the phase space. Consider a system described by a vector X(t) 
in configuration space. At time t w the system is cloned into a new system described by the vector Y(t). The two 
copies are let evolve with different realizations of thermal noise and the overlap Qt m (£) = X(t w + t) ■ Y(t w + t) is 
recorded as a function of time t. For coarsening-like systems (called type I systems inE3) the overlap converges to a 
finite value Qoo =_Limt^oo Qt m (t) f° r anv value of t w while for glassy systems with structural glasses behavior (called 
type II systems inEil) that limit gives, for all possible values of t w , the lowest possible value of X(t w + t) ■ Y(t w + t). 
Since in general it is possible to define IS also for coarsening systems, this is an indication that in these systems 
the non-equilibrium dynamics goes through configurations that are unrelated with the IS. Most probably the relevant 
configurations for coarsening are those on the border among IS-basins, i.e., those configurations which are not mapped 
into any IS. In this case the IS and the SW configurational entropy can still be defined but are obviously of little use 
for understanding non-equilibrium dynamics, as shown by the results reported in the following Sections. 

Before addressing the SW approach to non-equilibrium dynamics of constrained kinetic models we note that the 
definition of the SW configurational entropy may not be completely free of ambiguities, especially for systems with 
discrete states. Indeed the SW mapping assumes that once the energy and the dynamics are given then the mapping 
between configurations and IS is uniquely defined. This means that regardless of when a given configuration appaaxs 
in the dynamical evolution it will always be mapped to the same IS. In most of the recent papers on this subjectE3~E£l 
this was the case, but for discrete models with discretized values for the energies, as the ones studied here, there 
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may be problems because there could be many directions in phase space where the energy decreases by the same 
amount. In this case some "decision" must be taken, e.g., one could employ a random choice among the possible 
directions. The IS and the corresponding SW configurational entropy can still be defined, but now they depend on the 
chosen strategy for dealing with equivalent directions. It can be shownEZI that this leads to a temperature dependence 
of s c , so that the form of s c depends on T. Only when some additional requirements are fulfilled the temperature 
dependence disappears. For example when all possible strategies lead to the same sampling rate for the IS relevant 
for the dynamics. We anticipate that this requirement is not fulfilled by the constrained kinetic models nor by the 
BG model studied here, and in some cases we find different curves for s c for different temperatures. Nevertheless 
the violation is not too strong since we basically find only two different curves depending on the temperature range. 
This reflects the trivial fact that the constraint is more or less effective depending on the temperature because more 
we lower the temperature more the system orders. Another interesting feature of these models is that, despite the 
fact that the properties of these dynamical models are rather different, all of them have the same SW configurational 
entropy s c . This casts doubts on the relevance of the IS analysis for the non-equilibrium dynamics of these models. 



B. The kinetically constrained Ising chain 

These models are defined by eqs. ([[]) and (|^). To compute s c (e) we thermalize the system at a finite temperature 
T . This can be achieved either by running the dynamics for a sufficient long time or by starting from equilibrium 
configurations whose distribution is given by 

PeM) = (9) 

1 + exp(/j) 

For each thermalized configuration the corresponding IS is computed via the minimization process given by the zero- 
temperature dynamics. Repeating the procedure for several initial configurations the IS probability distribution (^) 
can be evaluated. 

Each IS is a fix-point of the dynamics, and therefore we can estimate s c from the number of fix-points. This is 
easily evaluated denoting with A the one-bit sequence I and B the two-bits sequence 01 since all fix-points are given 
by all possible arbitrary sequences of A's and B's (for instance, the sequence ABBBAABB). If Na,Nb stand for 
the number of A's and -B's in the sequence, with Na + 2Nb = N, where N is the length of the chain, then 

_ (n a + n b )\ 

Nfix - N A \N B \ ■ (10) 
From this expression, and noting that e = — (N A + Ng)/N, we have 

8 c (e) = l0g( j^ x) = -e log(-e) - (1 + e) log(l + e) + (1 + 2e) log(-l - 2e). (11) 

The above formula assumes that all IS are counted with an equal a priori probability. This is what is called the 
method of unbiased guess in information theorytJ. Therefore the above expression is valid iff the dynamics samples 
(almost) all fix-points with equal probability, where with "almost" we mean the fix-points relevant for equilibrium 
dynamics at temperature T. 

A better estimate of s c comes from the analysis of the zero-temperature dynamics. As shown in Appendix A the 
zero-temperature dynamics can be solved exactly for any value of a, and from this Pjs(e, T) can be evaluated: 

p ^vmmM J w^y (12) 

The details of the calculation are reported in the Appendices together with the expressions of (Co(oo)) c and (e/s), 
eqs. ( |6l| ) and (H|). Note that this result, as well as eq. ((l|), does not depend on the value of a implying that 
all models, and in particular the SCIC and ACIC, have the same Pjs(e,T) and hence the same SW configurational 
entropy. 

In Fig. [l] we compare the numerically evaluated Pjs(e,T) for the SCIC model with N = 64 and different temper- 
atures with the analytical prediction eq.(^). The agreement is quite good at low temperature but decreases with 
increasing temperatures where the variance of the Gaussian is slightly larger than in simulations. 
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FIG. 1. Probability histograms in the SCIC model with 64 spins at different temperatures compared with the analytical 
prediction (|l^) (full curve) and a fit to a Gaussian (dashed curve). 

In Fig. H we report the average IS energy as a function of T and compare it with the prediction from the analytical 
calculation of Appendix B. There are two possible ways to compute e(T): the first is from eq. ([59|), the second rests 
on integrating the variance 



e(T) 



(Co 2 M) c 

J^2 



dT, 



(13) 



where (Cq(oo)) c is given by (|6l|). We also report the result from the fix-point approximation (|iT|). Fig. § clearly shows 
that the different approximations depart each other at a temperature T ~ 0.6. Above this temperature the energy 
dependence of the IS free energy / in eq. (^) cannot be neglected anymore, showing that the fix-point approximation 
which neglects thermal fluctuations inside the IS-basins is inappropriate. On the other hand the direct calculation 
from the zero-temperature dynamics turns out to be very good. 
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FIG. 2. IS energies as function of temperature obtained integrating (pj|), using ( |59| ) and the variance of the IS-energy 
distribution given by (|l3|). 



Finally we consider s c . In Fig. (0) we show the results obtained for the SCIC model with N — 64 and different 
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temperatures. The SW configurational entropy s c is obtained from the numerical Pis(T,e) aaiSH 

s c (e) = 0e + log(Pi S (T, e)/N) + const. (14) 

For each temperature the constant has been fixed by collapsing different data onto the single curve. As a comparison 
we also show the theoretical predictions from equations (O) and [See Appendix B] 



. as) 



As shown in the Appendix B, both coincide asymptotically close to the ground state energy e = — 1. The collapse is 



excellent showing that the approximation (|TTj) and the low- temperature behavior (15) asymptotically coincide in the 
limit T — > 0. We note that there is a range of energies where data from T < 0.6 collapse on one curve while data for 
higher temperature collapse on a different curve. As discussed above, this residual temperature dependence follows 
from the presence of many equivalent directions for energy minimization. 

We have checked that Pjs(T, e) is independent of a by repeating the analysis for the ACIC model and for different 
values of a. In all cases we always find the same results. 




e 

FIG. 3. SW configurational entropy in the SCIC for N = 64 spins at different temperatures compared with the analytical 
prediction (|l^) (upper curve) and the fix-point estimate ([H]) (lower curve). 



C. The BG model 



In this case we cannot exactly solve the zero-temperature dynamics of the model and compute the SW configura- 
tional entropy. Nevertheless, we can approximate s c (e) of the BG model by counting the number of ways in which 
two or more particles can be distributed in a set of N occupied boxes separated by empty boxes. This yields two 
contributions: the first comes from all possible ways of distributing the M occupied boxes in a chain of N boxes, 
with the additional condition that each occupied box is surrounded by an empty box. This is again given by eq. jl0| ) 
assuming that a = 1 for occupied boxes and a = for empty ones. The energy (|5|) is given by E — —{N — M) and 
therefore this contribution reads, 

( e ) = -(1 + e) log(l + e) - elog(-e) + (1 + 2e) log(-l - 2e). (16) 

The second contribution follows from considering all different ways of distributing the particles among the M 
occupied boxes with the constraint that each occupied box contains at least two particles: 

1 N 
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where the n r \ terms arise from the distinguishability of particles. Introducing the integral representation for the 
delta function, 



1 f°° 

S(x) = — / exp(iAa;) dX, (18) 

2lT J-oo 

we find an expression for ( [l7| ) in terms of the fugacity y. In the N — > oo limit this can be evaluated by the saddle 
point method yielding 

logf N second ) 

s secon d[e) = = lQg{y) + {1 + e) log(exp(y) - tf - 1), (19) 

where y satisfies the saddle-point condition, 

c - l| cxpfa)-l~y _ 
y(exp(y) - 1) 

The full entropy is given by 

s c (e) = s f c irst (e) + s s c econd {e). (21) 



Note that for this model the configurational contribution may be negative because particles are distinguishable. 

We have computed s c (e) numerically following a procedure similar to that described for the constrained kinetic 
models. The results are shown in Fig. ^ for two different sizes N = 100, 500 and temperatures ranging from T = 0.1 
up to T — 1. Similarly to what found for the constrained kinetic models, the data nicely collapse onto a single curve 
although it does not exactly coincide with the number of fix points. In this model the presence of different equivalent 
directions to decrease the energy does not influence s c . This is most probably due to the global character of the 
constraint. 

Comparing Figs. [5] and ^ we see that the agreement now is worse. We attribute this to the presence of entropic 
barriers which follows from all possible arrangements of particles inside the boxes. All arrangements leave the energy 
unchanged, but their number strongly depends on the number of empty boxes, leading to a stronger energy dependence 
of the IS free energy for this model. This effect was not present in the kinetically constrained Ising chain. 




e 

FIG. 4. SW configurational entropy in the BG model for N = 100, 500 boxes at different temperatures 
T = 1.0, 0.5, 0.4, 0.3, 0.2, 0.15, 0.1 compared with the fix-point estimate @ (full line). 

The conclusion that can be drawn from this Section is that for these models a description of the glassy behavior in 
terms of a complex energy landscape is not relevant. Even though the SW configurational entropy for the constrained 
Ising chain is a non trivial quantity, it does not distinguish the SCIC model from the ACIC model. 
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IV. ANALYSIS OF COARSENING BEHAVIOR 



In the previous Section we have seen that the IS approach yields identical results for models which are known to 
have a completely different dynamical behavior, namely the SCIC and ACIC models. The purpose of this Section is 
to evidentiate these difference making connection with results already known in the literature and studying new ones 
to gain some insight using the tools from disordered systems. .-Coarsening appears when domains of a given phase 
grow in time slowly enough for the system to be off-equilibriumEa. In the simplest case, dynamics is characterized by 
a unique length scale associated with the typical size of the growing domains. All models discussed in this paper can 
be described in terms of coarsening in the sense that it is possible to define a length scale which identifies the distance 
from to equilibrium. For the kinetically constrained Ising chain this length is the typical size of the +1 domain while 
in the BG model it is the typical length of sequences of empty boxes. 

In the simplest cases this length suffices to characterize the off-equilibrium behavior. For instance, for coarsening 
in (ordered or disordered) ferromagnets the off-equilibrium behavior is fully characterized in terms of a single length 
scale L(t) in the sense that the two-times dependence of correlation and response functions directly enters through 
the value of this length scale. In the aging regime, where both times t, s are large we haveEj 

(22) 



~ 777v^(tS)> t>s > ( 23 ) 



with a > 1 a positive exponent which depends on the model under consideration. 

Correlations are easy to measure, while responses require the introduction of an external perturbation. This 
perturbation must couple with the variables of the system and must be small enough to ensure a linear response 
regime. A useful quantity is the integrated response function (hereafter denoted as IRF), pinch measures how much 
the system remembers the effects of a perturbation applied during a given interval of timeE£l, 

X(M) = / duR{t,u) (24) 

J S 

One of the salient features of coarsening phenomen is that aging effects in the integrated response function 

asymptotically vanish when the lower time s goes to infinity and the system does not have long-term memory. Suppose 
that the coarsening length L(t) grows like t 1 ^ with z a dynamical exponent. Using eqs.(^3|) and fl24| ) one finds that 
the aging part of the IRF behaves like, 

X ag (t,s)r.s 1 -^x(t/s) (25) 

which vanishes as s — > oo if a > z. 

An easy way to test these effects is by directly looking at fluctuation-dissipation plots (hereafter referred to as 
FDT)c3. In equilibrium R(t,u) = R(t —u) = f} ° g~ which substituted into ( pi|) yields, 

X (t-s)=0[C(Q)-C(t-s)) (26) 

and the plot of T\(t, s) in terms of C(t, s) is a straight line of slope — 1. _ . 
In the off-equilibrium regime expression (^6|) can be generalized by defining the fluctuation-dissipation ratioE^I 

TR(t,s) , , 

ds 

which measures how far the system is from equilibrium. In equilibrium X is equal to 1. In the off-equilibrium 
asymptotic long-time regime, i.e. in the aging regime, where there is no time translational invariance both X(t,s) 
and x(ty s ) are expected to be non-trivial functions of C . A quantitative estimate of A, can be obtained from the 
slope of the FDT plots: 

X(C) = - [ 9T X(M) 1 _ (28) 

. 9C{t,s) \ C {t,a)=G 
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For coarsening models the aging part of the IRF asymptotically vanishes and x(i, s) is expected to saturate to a 
finite value, called the field-cooled value in the context of spin glasses, and stays constant while the correlation still 
decreases before saturating. 

The simplest way to compute the IRF in the present class of models is to apply a perturbation which does not 
couple with the absorbing state, i.e., the ground state. We start from a random initial configuration and after a 
waiting time s = t w apply a random staggered field hi — ho(t), where ho(t) is the intensity of the field and e; = ±1 
are independent random quenched variables of zero mean. This method has the advantage that the perturbation term 
in the Hamiltonian does not directly couple to the coarsening length, and has been used by BarratCJ to investigate 
coarsening in finite-dimensional Ising models. 

The staggered magnetic field couples to the spin variables Ui in kinetically constrained models and to the equivalent 
variables ct; = 5 nu o in the BG model so that the perturbation in the Hamiltonian reads 

N 

6H{t) = -'%2e i k (t)<Ti . (29) 
»=i 

In the case of ho(t) = ho 0(t — t w ) the integrated response function can be obtained measuring the random staggered 
magnetization after switching the field at time s — t w as, 

1 N 

^ i — 1 

The original 0, 1 variables have some disadvantages, for example the correlation at equal time is not 1 but depends 
on temperature. For this reason we find more convenient to work with the new variables v = 2a — 1 which now take 
the values 1,-1. We then consider the disconnected correlation, 

1 N 

C(t + t w , t w ) = — Vj{t w ) Uj(t + t w ) (31) 

i=l 



and the staggered magnetization (p0[), 

1 N 

M stag (t + t w ,t w ) — — 2J £j Vi(t + t w ) = 2hox(t + t w ,t w ) (32) 

8=1 

Now the equal times disconnected correlation function is equal to 1 so that in the FDT plots, where the integrated 
response function is plotted versus the disconnected correlation function, all curves start from C — l,M stag — for 
t = t w . This makes easier to compare results from different values of t w . 

As discussed before in Section III a different way to distinguish coarsening dynamics from other more complex 
behaviors is to measure the overlap Q{t) between two replicas which start from the same configuration at time t w and 
evolve with different realization of thermal noises. Here we consider the connected, normalized overlap function, 



m _ j EL 4 1] fa + t) og (t w + t) EL ^ (jjo+mh Xli ^ fa + «)) ,™ 



where ct^ 1,2 -* = 0, 1 refer to the two replicas, and the connected, normalized, correlation 

<~> C {1 -t l W ,l W ) — — 7TT " — 7TT ^04J 

For type I systems — limt^oo Q% {t) is finite. For type II systems this quantity converges in the t — > oo limit 
to the lowest accessible value, i.e. vanishes in the t — > oo limit. In equilibrium both ( [33] ) andJ0) depend only on 
t and the following relation is valid Q% (t) = C c (2t). There are few numerical studies of Q^EjO. For the models 
considered here the results from this analysis are, however, not so strong, the reason probably being that coarsening 
occurs in a disordered (i.e. paramagnetic) phase. 

In the rest of this Section we investigate in detail coarsening length scales, correlations and responses for the 
non-equilibrium dynamics of models described in Section IL_iMc_note that although many results on timescales and 
coarsening length scales have been obtained in the literaturetlltjEZI almost nothing is known about the aging behavior 
in this type of models (partial results are shown ir£3 for the SCIC). 
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A. The SCIC model 



It has been shown ini that the SCIC model has an activated timescale t± — exp(/3) characterized by an exponential 
decay of the staggered energy. For times smaller than t\ there are no aging effects and only for times larger than t\ 
non-equilibrium behavior with non-exponential relaxation and aging appear. From the decay of correlation functions 
a second activated timescale r corr > t\ can be defined. To this end we have computed the equilibrium connected 
correlation function 

r m nzZi ^(Q) *m - Etx <m EL ™ 

N 22i=i cri(0) - (Ei=i ^(0)) 2 
which is well described for all temperatures by the following functional form 

C &) = ( i ,\, ) exp(-at fo ), (36) 



l + t/n 

where a, o, b are three fit parameters and the activated timescale T\ is introduced in the fitting function as an effective 
microscopic time. The results for C c (t) are shown in Fig. |^ for different temperatures, the lines are theJaest fits with 
form (p6|). This fits arc in agreement with the asymptotic analytical predictions of Reiter and JackleE3 and Schulz 
and Trimpercfl but combiii6|ci-svith the exponential timescale T\ derived inB3. In particular the exponent a is close to 
the value 1/2 predicted irO'Ell for very low temperatures. 




From the fit we can estimate r corr as 

oo 

C c (t) dt . (37) 

The results are shown in Fig. |6j. The correlation time r corr follows the Arrhenius type law r corr ~ exp(2/J). This 
functional dependence of correlation time from temperature can be understood from the following phcnomcnological 
argument, based on defects annihilation in the SCIC. A defect separated by magnetized domains can disappear by 
anchoring defects along the chain. The typical time to anchor a defect is exp(/3) while the length of the magnetized 
domains in equilibrium is of order exp(/3). Because defects can be anchored starting form the right or from the left 
of a magnetized domain, the typical time to annihilate that domain is the sum of independent processes yielding 
~exp(2/3). 
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(3 

FIG. 6. Correlation time in the SCIC computed using (^) and fitted with an Arrhenius behavior. The best fit gives 
r corr = 1.43exp(1.93/3) 



To investigate coarsening in the SCIC model we have measured the growth of the average domain length 



where 



N i+l-1 
i=l j=i 



(38) 



(39) 



is proportional to the probability of having a domain of spins 1 of length / at time t. For t — ► oo Pt converges towards 
the equilibrium length probability distribution 



^(0 



exp(-/3) 



[l + exp(-/3)]' 

and the average length saturates to the equilibrium value 

d oq = 1 + exp(/3) 



(40) 



(41) 
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FIG. 7. Average domain length and magnetization in the SCIC model. The average length grows diffusively like t 

In Fig. |?] we present the average domain length as a function of time starting from a random initial condition. From 
the figure it follows that after 10 6 MCS d(t) is still well below the equilibrium value d cq for temperatures as high as 
T = 0.2 indicating that the. systems has not yet equilibrated. In agreement wither a power law fit leads to d(t) ~ t 1 / 2 
characteristic of diffusionc3. In the lower part of Fig. |t] we show the relaxation of the energy as a function of time. 
These results combine the zero temperature exponential decay to the threshold energyBZIcfl 1 /e with the slower decay 
towards equilibrium. Note that since the average domain length grows like i 1 / 2 while the equilibrium value for large 
(3 is d eq ~ exp(/3) we get for the correlation time r CO rr ~ exp(2/3), as expected. 

Further informations on the non-equilibrium behavior can be obtained from the analysis of the response to a 
staggered magnetic field as described in eqs. (p8|), d|l]) and (^2|). In Figs. || and § we show M stag (t + t w ,t w ) [see 
eq. (p2[)] and the correlation function C(t + t Wl t w ) for temperature T = 0.3 and 0.11 and different waiting times t w . 
The strength of the staggered field is ho = 0.1 while the system size is N — 10 5 . The horizontal lines indicate the 
equilibrium values 



C e 



1 N 



1 - cxp(-/3) 
1 + exp(/3) 



(42) 



and 



N 



/?exp(/3) 
[l+ex P (/3)] 2 



0{hl) 



(43) 
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FIG. 8. Correlations and zero-field cooled magnetization in a staggered field in the SCIC model for N = 10^, T = 0.3 and 
different waiting times t w = 10, 100, 1000, 10000. The horizontal lines indicate the equilibrium values ([I2]) and (fl3|). 
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FIG. 9. Correlations and zero-field cooled magnetization in a staggered field in the SCIC model for TV = 10 J , T = 0.11 and 
different waiting times t w = 10, 100, 1000, 10000. The horizontal lines indicate the equilibrium values (|l3"|). 

The corresponding FDT plots are shown in Figs. [h] and |ll|, respectively. 
Concerning Figs. H |, |g and [fl] we note the following: 

1. Aging in the correlation function appears for waiting times larger than the critical time n = exp (/3) and 
survives even for times larger than the correlation time r corr . This can be seen from both correlation function 
and staggered magnetization. The equilibration time r cq to reach equilibrium is larger than both the activated 
time T\ and the correlation time r corr . If we define the equilibration time r cq as the time needed for the average 
domain length to reach the equilibrium value, then plotting the data of Fig. |?] as function of Tlog(t/ro) we find 



that T r , 



Tl 



i.e. t. 



cq 



exp (3/3). 



2. The staggered magnetization has a hump, in correspondence of which the correlation function presents a broad 
minimum, as function of t. For t w < t\ the hump maximum takes the largest value and decreases with t w as 
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soon as t w > n and eventually disappears for t w — » oo. This effect is a direct manifestation of the two critical 
timescales present in the SCIC modelEZI. 

3. The existence of different activated relaxation times results in rather peculiar FDT plots, see Figs, [h] and [n|. 
For t w < t\ , Fig. [HJ C, x an d X do not show any dependence on t w , nevertheless X is a non-trivial function of 
C corresponding to non-equilibrmpi behavior without aging. A similar shape is found for the one-dimensional 
fsing model at low temperaturesej. For t w > n, Fig. [l(j, there are aging effects and X shows the typical two 
slope pattern. However, the existence of a second typical timescale results in a second downwards bending of 
the IRF and X as function of C has a three slope shape. 
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FIG. 11. FDT plots in the SCIC for N = 10°, T = 0.11 and different waiting times t w = 10, 100, 1000, 10000. The straight 
line is the FDT relation @ 

We conclude the analysis of the SCIC model by discussing the results for the cloning experiment. In Fig. |l2| we 
show the overlap Q\ {€) [eq.(fS3|)] as a function of the connected correlation C c (t + tw,t w ) [eq.(p3)] for temperature 
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T = 0.11. The results show that, for any finite t w , Q goes to zero quite rapidly. If one compares Q c t (t) with 
C c (2t + t w ,t w ) for different values of t w one finds that both decrease exponentially fast with time and that Q£ (t) 
is smaller but very close to C c (2t + t w ,t w ). This implies that the two trajectories depart form each other quite fast. 
The data in Fig. [13 collapse quite nicely onto the parabola Q ~ C 2 in agreement with the exponential decay. 
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FIG. 12. Q c tw {t) versus C c (t + t w ,t w ) in the SCIC at T = 0.11 



In conclusion the SCIC model is a coarsening model with three activated timescales T\ = exp(/3), r corr = exp(2/3) 
and r oq = exp(3/3). In this scenario we do not expect the simple scaling form (|2^) in terms of a single length scale 
to be verified. Indeed, although a growing len gth scale can be identified a single length scale does not describe the 
whole time regime. This is clearly seen in Fig. hj for T = 0.15 where aging starts after t\ ~ 1000. The scaling form 
( p2| ) constructed from data from Fig. |t] is shown in Fig. [l4|. The scaling is obviously rather poor. The behavior 
of this model is essentially diffusive as emerges from the behavior of the overlap Q c t it) similar to the behavior of a 
ferromagnet with non conserved order parameter, with the difference that the SCIC model has no phase transition 
and coarsening takes place in a disordered phase. 
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FIG. 13. Connected correlations C c (t m ,t m + i) in the SCIC for iV = 10 5 , T — 0.15 and different waiting times 
t-w = 10, 100, 1000, 10000. Aging is present for times larger than n = exp(/3) 
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B. The ACIC model 



Coarsening in this model has been extensively studied tw-,saveral authors, finding that the correlation time has 
a super-Arrhenius behavior and grows like exp(/3 2 / log(2))c3'E3, much faster than the typical Arrhcnius behavior 
~ exp(2/3) found for the SCIC model. In the SCIC model domains can always grow if they can annihilate defects by 
building intermediate defects to the right or to the left of that defect. In the ACIC model, on the contrary, defects can 
disappear only by anchoring intermediate defects in the middle of the domain from one side. This strongly enhances 
the correlation time. On the other hand, the absence of a critical time like n and the coincidence of the correlation 
time with the equilibration time makes the dynamics of this model simpler than that of the SCIC model. 

In Fig. |l^ we show the average domain length defined by eq. (|3|) and the energy as a function of time when 
startiHgi from a random initial configuration for different temperatures. The Tlog(i) scaling predicted by Sollich and 
Evansc3 is very well satisfied. Note the presence of plateaus in both the average domain length and the energy for the 
same range of time. These correspond to time intervals where domains coalesce and the global eaergy stays constant 
because the number of anchoring spins is much smaller than the length of the coalescing domainsE3. Since the average 
domain length grows like log(d) ~ ATlog(i), i.e., d ~ t XT , and the equilibrium domain length is given by d cq ~ exp(/3) 
for low temperature, a phenomenological argument yield&-for the correlation time r corr ~ exp(/3 2 /A) with A = log(2) 
in agreement with the expectations of Mauch and JacklcEj. 
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FIG. 15. Average domain length and magnetization in the ACIC at T = 0.15, 0.20. The average length grows like t 

Figures [l6] and |l^ show the non-equilibrium disconnected correlation function ( |3l| ) and the staggered magnetization 
( |32| ) from infinite temperature initial conditions as a function of time for two temperatures T = 0.2 and 0.4 and different 
values of t w . The strength of the field is ho = 0.1. The dashed horizontal lines are the equilibrium values (42) and 
@. The FDT plots are shown in Figs. [l| and |l|. 




t 

FIG. 16. Correlations and zero-field cooled magnetization in a staggered field in the ACIC model for N = 10 5 , T = 0.4 and 
different waiting times t w = 10, 100, 1000, 10000. The horizontal line indicates the equilibrium value (k3h 
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FIG. 17. Correlations and zero-field cooled magnetization in a staggered field in the ACIC model for N = 10 5 , T = 0.2 and 
different waiting times t w = 10, 100, 1000, 10000. The horizontal line indicates the equilibrium value (0) 



Looking at this set of figures we note the following points: 

1. The staggered magnetization does not posses a hump, and the correlation functions a broad minimum, as in 
the SCIC model. Here, on the contrary, both are monotonic functions of time, a behavior commonly found in 
models where there is no critical time (like T\ in the SCIC) associated with a microscopic fast process. On the 
timescale of r corr both quantities relax to the equilibrium values. We note that r corr ~ 10 4 for T = 0.4 and 
r corr ~ 10 15 for T = 0.2. 



2. A simple look at Figs. |l^ and 17 reveals that aging is present for all timescales in both the correlation function 
and staggered magnetization. Aging in M stag is noticeable for all values of t w suggesting that aging in the IRF 
disappears rather slowly with t w . Keeping in mind the coarsening nature of this model this suggests that a ~ z 
in eq. (|2T 



3. From Figs. and [l?] it is difficult to verify the coarsening nature of the dynamics in this model. A simple 
check can be done with the help of the FDT plots shown in Figs, [l^ and 19 for temperatures T = 0.4 and 
T = 0.2. Interestingly for waiting times comparable to the correlation time, so that the system is not far from 
equilibrium, the fluctuation-dissipation ratio X rapidly converges to 1, see Fig. [l8| At low temperatures, Fig. 
[To| , t w <C T corr and the fluctuation-dissipation ratio is very small, X ~ 0.1, and roughly independent of t w , a 
scenario typical of coarsening modelsEj. 
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In the ACIC model there is only one characteristic divergent timescale, namely r corr , and thus we expect that the 



scaling behavior (|22j) should be satisfii 
grow in time like a power law [see RefJ 



Note that, contrarily to the SCIC, the average domain length d(t) does not 



and Fig. h5| . Consequently, in the aging re 



form t/t w , see Fig. 2^, but a more complicated function depending on shape of d(t)t 
is plotted as a function of d(t w + t)/d(t w ). The scaling is quite good. 



sime the scaling will not be of the 
il, see Fig. f2l| where C c (t w ,t w +t) 
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FIG. 20. Connected correlations C c (t Wy t w + i) in the ACIC for N = 10 5 , T = 0.20 and different values of t w plotted versus t 
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FIG. 21. The same data of Fig. ^ plotted versus d(t w + t)/d(t w ) 

We conclude our analysis of the ACIC model with the discussion of Q% (t) and C c (t + t w ,t w ) shown in Figs. and 
p3| . Again Q% m (t) decays to zero for any t w but slower than in the SCIC model, see Fig. [l^. If (t) is compared 
with C c (2t + t w ,t w ) for different values of t w one finds that again Q% (t) is smaller, but very close to, C c (2t + t w , t w ). 
However, contrarily to the SCIC case, now during the initial regime when t w is small both Q° (t) and C c (2t + t w , t w ) 
show a plateau, more pronounced for Q, for times t ~ 10 3 . The time range between t w — 10 2 and t w = 10 3 corresponds 
[see Fig. [l5| to the regime where d(t) is growing very fast. This means that during this time interval domains grow, 
since C slowly decays, but Q remains almost constant because the two replicas follow the same narrow path in phase 
space. This effect is consequence of the way domains grow in this model where the anchoring of spins proceeds one by 
one in a given direction, different from the diffusive mechanism in the SCIC model. For t w larger than 10 3 this effect 
would be observed at the next timescale, between 10 5 and 10 6 , where new domains would have grown again [see Fig. 

Moreover this would also lead to a new plateau for Q and C for waiting times t w = 10 and 10 2 for values of the 
correlation of order 0.3 (which we have not reached in the simulations). 
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FIG. 23. Q c tw (t) and C c {2t + t w , t w ) in the ACIC at T = 0.2 

We can summarize the results of this subsection by saying that the non-equilibrium behavior of the ACIC model 
resembles coarsening in a simple ferromagnet although the growing length scale grows slower, like t Tlo ^ 2 \ The 
ansatz (^2|) for the scaling behavior is well verified and the FDT plots show some similarities with the physics of 
undercooled liquids although no connection between the slope of the FDT behavior and the SW configurational 
entropy is possible. Note that we obtain X ~ when temperature is lowered. Furthermore, the behavior of Q c t (t) 
reveals that the relaxational dynamics proceeds by evolution in narrow channels in the time intervals where domains 
grow similar to what happens for type I models. On the contrary, it is not possible to establish the glassy scenario 
from the value of limt_ >c>0 Q\ (t) a quantity which we expect to work only for models with broken ergodicity. 
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C. The BG model 



In the ID BG model, at difference with the kinetic constraint models discussed above, there are no dynamical 
constraints and coarsening follows from the slow growth of the number of empty boxes induced by entropic barriers. 
If we denote with a = 0, 1 empty and occupied box, respectively, then we can consider the same quantities discussed 
for the kinetic constraint models, i.e., domain length probability distribution, average length, correlation and mag- 
netization. All these quantities converge for large times to their equilibrium values derived from the equilibrium 
probability distribution^, 

^(n) = lf(^, n )= Z "' X , eXP( f^ 0) (44) 
N ' n\ exp (z) 

r—l - 1 

which gives the probability that at equilibrium a box contains n particles. Normalization of the distribution ( |4^ ) 
corresponds to the conservation of the total number of particles, and reads 

exp(/3) - 1 = (z - 1) exp(z). (45) 

By using P oq the equilibrium values of correlations, response and domain length can be computed. For example the 
average domain length of empty boxes at equilibrium is given by, 

d - " idyo) (46) 

which diverges for f3 — > oo when all particles fill a single box. From the definition it follows that P e q(0) is equal to 
minus the equilibrium energy (j|) which, for small T, goes as E/N = — 1 + T + 0{T 2 ). We then have d cq ~ for 
/3->oo. ' 

Similarly for the correlation function we have, 



1 N 

C eq = ^^(^) 2 = [2P eq (0)-l] 2 (47) 



i=l 



where v — 2a — 1 . 

To compute the staggered magnetization we have to use the equilibrium probability distribution with the extra 
term (^9|) added to the energy. A straightforward calculation leads to 

M cq _lf , v _ efl(l-e-') smhQ3V) 

M stag - N 2_. W - z{e z_ 1 + e C Osh(/% )) ' 



where z is now solution of 



ze 



i i 



1 + e /3(l+M e z _ 1 + e 0(l-h o ) 



(49) 



Note that M s c t q nj is linear in h for small values of h . 

Analyzing in details the dynamics we can distinguish two different decay processes: the first one entropically 
activated and the second one energetically activated. When the system is quenched from high to low T particles 
"evaporate" from some boxes and accumulate in others. After some time a situation is reached where boxes with 
more then one particle are separated by empty boxes and a very small number of single occupied boxes (defects). For 
T = these defects disappear and the energy does not relax to equilibrium. On the contrary for T small but finite 
the number of these defects may be large enough, its number scaling as TN, to serve as nucleation paths between two 
nearby multiple occupied boxes which, by the usual entropic mechanism, eventually accumulate onto a single box. 
The timescale T\ for this process is activated since defects must be created, but it is smaller than exp(/3). 

A second, energetically activated, process appears for defects to be anchored between multiple occupied boxes so 
that they coalesce in a single multiple occupied box. At low temperature and close to equilibrium the typical number 
of particles per occupied box is of order /3 while the distance that particles must cover by diffusion from one box to a 
contiguous one is of order d cq ~ (3. Combining these two behaviors we obtain for the equilibration time r cq ~ /3exp(/3). 



The interplay between these two mechanisms can clearly be seen in Fig. 24 where we report growth of the average 



domain length, from an initial random configuration, and (minus) the average number of empty boxes (energy) as 
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a function of time for different temperatures. The data are plotted as function of Tlog(t) so that the equilibration 



time, where the average length reaches the equilibrium value, is Tlog(r c , 



1 + 8 with 8 = -Tlog(T) leading to 



This equilibration timescale r eq is shorter than that of the kinetically constrained Ising chain but 



T oq ~ /3exp(/3). 

longer than the mean-field case r cq ~ "y P ' B. From the figure we also see that both the average domain length and 
the energy display a plateau at short times corresponding to the T — behavior. The departure from this initial 
plateau occurs for times shorter than the activated characteristic time exp(/3) and is driven by the entropic mechanism 
described above. A good collapse of the departing time is obtained for n = exp(/3)//3. For (3 — > oo the characteristic 
times Ti and r eq become well separated since r cq /ri ~ 1 . 




Tlog 10 t 



FIG. 24. Average domain length and energy in the BG model. 

Figures |2?] and [2^ show the correlation function, staggered magnetization and FDT plots of the ID BG 

model for the two temperatures T — 0.1 and 0.09 and field strength ho = 0.1. 




t 

FIG. 25. Correlations and staggered magnetization in the BG for iV = 10 4 , T = 0.10 and different values of t n 
From the analysis of the figures the following conclusions can be drawn. 
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1. After quenching from infinite to low temperature a fast evaporation of occupied boxes occurs after which only 
a finite fraction of them, approximately 30%, survives. Each occupied box contains in average about 3 — 4 
particles. This process is clearly seen in Fig. ^4|. The timescale t\ ~ exp(/3)//3 has a role similar to the timescale 
Ti found in the SCIC model. The aging effects are absent for t w < t\ but X < 1 [see Figs^27| and ^8). The 
dynamics is again diffusive and similar to the one found for the one-dimensional Ising modelEl 

2. For waiting times t w > t\ the dynamics slows down due to energy barriers. To empty a box a large number of 
particles must be transferred, a process which is cooperative and involves all the particles in the box. The typical 
time of this cooperative process is r oq ~ /3exp(/3). For waiting times t\ < t w < r eq the system shows strong 
non-equilibrium effects with a downwards bending of the IRF as a function C similar to what seen in the SCIC 
model. The origin of this effect is, however, different and follows from the asymmetric response to the staggered 
field of occupied and empty boxes. Since the field is coupled to empty boxes, the typical time to empty a box is 
larger than that to occupy an empty one. In other words, when quenching from high (or infinite) temperature 
boxes are occupied fast and its number converges relatively fast towards the equilibrium value. However, due 
to the staggered field, the distance among them is far from the equilibrium value and occupied boxes must be 
rearranged, which is a very slow process. Consequently correlation functions and staggered magnetizations show 
peculiar humps corresponding to the fast and slow responses. 




FIG. 26. The same as Fig. for T = 0.09. 
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In conclusion we can say that the strong entropic effects which follows from the large occupation numbers of some 
boxes at low temperatures, imply that the non-equilibrium behavior of the ID BG model cannot be well described 
only in terms of a single length scale d(t). The reason is that although d(t) tells how far the system is from equilibrium 
it does not contains enough informations to efficiently describe the effects associated with the entropic barriers and 
the ansatz ( p2| ) docs not hold. 



V. CONCLUSIONS 



In this paper we have studied the dynamics of constrained ID Ising models. In particular we have focused on 
a family of constrained kinetic models which interpolate between the symmetrically constrained Ising chain (SCIC) 
introduced by Fredeickson and Andersenlij and the asymmetrically constrained Ising chain (ACIC) introduced by 
Eisinger and Jacklccil. For comparison we have also studied the ID Backgammon modeled where dynamics is slowed 
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down by the global constraint imposed by the conservation of the particle number. Although these models reproduce 
some generic features of undercooled liquids there are, however, important differences. 

First of all since the class of kinetic constrained models defined in Section II, [see eqs. ([!]) and (0)] have the same 
Stillinger- Weber configurational entropy [see Section III] but rather different non-equilibrium behaviors [see Section 
IV], we conclude that the IS approach is not appropriate for these models. The reason of this failure can be easily 
understood. Dynamics in coarsening systems is usually described in terms of a growing length scale which measures 
how close the system is to equilibrium. In this scenario configurations at different times are obviously overlapping 
since domains at the early times are contained in the larger domains at later times. Consequently the (slow) approach 
to equilibrium can be described within a geometrical picture in terms of the value of the average domain length. 

In undercooled liquids such a length scale probably does not exist and no coarsening takes place in the metastable 
region. Here the slowing dewn of the dynamics follows from an activated dynamics in a complex energy landscape 
composed by many valleysc3. The SW idea of mapping run-time configurations onto local minima of energy surface 
(IS) is to provide a statistical description of the valleys. The fact that these valleys are uncorrelated is at the basis 
of the potentiality of the approach, a result .shared by disordered mean-field models for the glass transition such as 
p-spin modelsc3 or the random energy modelES. Obviously coarsening plays a role for crystallization processes butjw: 
know that the anomalies found in the undercooled regime are also found in disordered models with a crystal stateBfl, 
so that the role of the crystal state in the undercooled dynamics can be ignored. 

We have critically discussed the SW configurational entropy and its meaning. A conclusion which_.is_difficult to 
escape from is that any sort of configurational entropy, for example adapted from mean-field theoryoEZI, will have 
meaning only from a dynamical point of view. Efforts in this direction arelhc Stillinger and Weber approach itself, 
the work of NieuwenhuiacrO, the mean-field scenario by Franz and VirasoscEU and the very recent approach proposed 
by Biroli and KurcharH. 

The description of non-equilibrium dynamics in terms of configurational entropy is valid, in general, if relaxation 
proceeds via activated jumps between uncorrelated configurations not described by any characteristic length scale. 
Another possible rephrasing of this conclusion, following the definitions inc3, is to say that models of type I cannot 
be described in terms of any mean-field like configurational entropy such as the SW approach. Inherent structures 
and their statistical treatment are only useful for models of type II, category to which structural glasses belong. But 
as we have shown here, one must be careful in using this classification for models without phase transition but with 
glassy dynamics when ergodicity is not broken like in ordered (or disordered) ferromagnets. A simple classification 
based on the infinite time limit of the Q tul (t) seems to be well posed fcp-|Coarsening ferromagnet-like systems, but 
for other slow dynamics it may be inappropriate, see some results in Ref.E-3. Still a careful examination of the whole 
dynamical behavior of Qt m (t), and not only its infinite time limit, may provide useful information to evidentiate how 
coarsening takes place. 

Complementary to this investigation we also made a detailed study of the coarsening behavior of the ID constrained 
kinetic models and the ID Backgammon model analyzing their different non-equilibrium behaviors. Such study 
confirms theoretical results already present in the literature, but also contains new ones concerning the behavior of the 
average domain length, in particular for the ID Backgammon model for which nothing is known, of the magnetization, 
responses and the violation of fluctuation-dissipation relation. Concerning correlation functions we have found that 
in all models studied here the t/t w scaling behavior characteristic of mean- field models is not observed. Howeuer, 
when there is a single typical time or length scale, as in the ACIC model, the_general scaling ( p2] ) and (^3])E3 is 
found. For the SCIC model we have found that the predicted critical timescale TiE3 marks the violation of the scaling 
behavior (|2^). This peculiar behavior leads to a hump and a broad minimum in the integrated response function and 
the correlation respectively, leading to a rather unusual fluctuation-dissipation plots. This scenario remains valid for 
values of a not too far from 1/2, the region of a close to 1 (or 0) is under investigation and will be reported elsewhere. 
A similar conclusion is valid for the ID Backgammon model where the presence of entropic barriers induces a second 
characteristic timescale. In this respect the model with the simpler non-equilibrium dynamics is the ACIC which 
presents a single timescale and the scaling behavior ( p2| ) is well satisfied. This results suggest that the FDT plots 
to distinguish one glassy scenario from another one must be used with caution. The simple models studied here 
show interesting FDT plots but no relevant informations can be obtained to interpret the off-equilibrium scenario. 
Quite probably the only cases where quantities such as configurational entropy or FDT ratio are interesting and 
meaningful are those where some kind of universality is expected. To this class belong structural glasses. For this 
case, the Stillinger- Weber decomposition provides a statistical and useful description of some, otherwise inapplicable, 
mean-field concepts. 

In the present research we have tried to clarify the limit of validity of the SW approach to the description of the 
glassy dynamics. The conclusions that can be drawn from the present study of kinetic constrained models and the 
Backgammon model can be generalized and we expect that the SW approach will fail to describe any coarsening 
model. In coarsening models, the information is always contained in a growing domain length which specifies the 
stage of evolution of microscopic configurations. In undercooled liquids where the simplest coarsening description 
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fails, the SW description is a good alternative which provides the natural link with mean-field ideas. 



ACKNOWLEDGMENTS 

We acknowledge P. Sollich for useful discussion. (F.R.) has been supported by the Spanish Government through 
project PB97-0971. MS is supported by the European Commission (contract ERBFMBICT983561). 

Note added: After completion of this work we have known of a workE2l by J. P. Garrahan and M. E. J. Newman on 
a three-spin interactions model on a triangular lattice who find results similar to those presented here. 

APPENDIX A: CLOSURE OF DYNAMICAL EQUATIONS 

For the constrained kinetic Ising chain the single-spin zero-temperature dynamics s: 

h = -Ti(aT i+ i + 6tj_i) (50) 
with a, b such that a + b = 1. For k > we get 

-jz(n ■ ■■ T i+k ) = -a{kn ■ ■ ■ n +k + t;_i • • -n +k ) (51) 

-b (k Ti ■ ■ ■ T i+k + Ti ■ ■ ■ T i+k+1 ) (52) 

FollowingEl we define the set of correlations, 

1 N 

Ck ^ = n I>W r '+i(*)- T '+fc(*) ( 53 ) 

It is simple to see that, for any a, b provided that a + b = 1, the C k satisfy the equations 

^l = -kC k (t)-C k+1 (t) (54) 
at 

which are identical to those of FA model computed i nE. The generating function 

OO 

closes the hierarchy yielding the result0, 



k\ 

k=Q 



G(x,t) =G ((l+a;)exp(-i)-l) (55) 
with the initial condition Ga(x) = G(x, t = 0). In the t — » oo limit ( p5| ) yields C k (oo) = Co (00)^0 where, 

00 /_-i\k 

k=0 



APPENDIX B: ANALYTICAL CALCULATION OF P IS (E,T) 

Here we give an analytical derivation of the SW configurational entropy for the constrained Ising chain from the 
zero-temperature dynamics described in Appendix A. If at t = the system is in thermal equilibrium at temperature 
T and the system is large enough the central limit theorem says that the C k (0) are Gaussian distributed variables. 
Consequently, Co (00) is a Gaussian distributed variable with the first two moments given by, 

^M) = E i ir^( )) ( 57 ) 

fc=0 
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00 ( i \k J rl 

(C7 2 (oc)) c = <C 2 (oc)} - (Co(co)) 2 = J2 ^Jji-(Ck(P)Ci(p))c (58) 



fell! 

fc=0 



where the subindex c stands for connected correlations. If m stands for the magnetization m = 1/[1 + exp(— /3)] 
then Cfc(O) = (1 — m) fc+1 for a thermalized initial condition at temperature T. This automatically yields, using ( |57| ) 
the average IS energy (e/s) = (Co(oo)) — 1 yielding, 

(e JS ) = (1 - m) exp(m - 1) - 1 . (59) 



The computation of the second moment (58) is also straightforward and requires the calculation of (CfcC/) c . A 
simple calculation gives, 

{C k C l ) c = ±-((k-l-l + -)(l-m) k+1 -(l + k + l + -)(l~m) k+l + 2 ) . (60) 
iv V m m / 

This expression is valid for k > I and k + I + 1 < JV. After some calculations we finally obtain, 

,„9/ . . If, . 4(1 — to) 2 , , (2 — m)(l — to) , , 

C oo c = - 1 - m 3 exp m - 1 - -i '- exp to - 1 + -/o 2(1 - m 

iv V " m m 

771 + 2 \ 

+ 2(l-m) 3 exp(2m-2) - ( — — )(1 - m) 2 exp(2m - 2) . (61) 

TO / 



where Iq(x) is the zeroth order modified Bessel function Iq(x) — X^fe^o 2 2 ^(fc!) 2 • This finally yields for the IS 
probability distribution (||): 

From eqs. ( ^9| , |6l|) we can directly obtain the configurational entropy in two different ways. One is obtained by exact 
integration of the IS energy as function of temperature, 

f T de(T) dT 

where e(T) is given by the expression (p9[). The other is obtained by integrating the fluctuations, 

f T (c 2 M) c 

S Jo T3 dT ( - 64 ' ) 

According to (||) both expressions should coincide when the term f(/3, e) in (^) is independent of the energy e. 
This is generally not true and expressions (|l5|), ( |64| ) are different. The interesting remark is that, in the limit T — > 
both expressions (|l5[), ( |64|) and the fix-point approximation ([ll|) coincide to first order in T. The reason is that the 
SW configurational entropy has full meaning in the limit where s c goes to zero and the IS-basins are very narrow 
containing few configurations. 
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